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Abstract 

We compare two definitions of multistable Levy motions. Such pro- 
cesses are extensions of classical Levy motion where the stability index 
is allowed to vary in time. We show that the two multistable Levy mo- 
tions have distinct properties: in particular, one is a pure-jump Markov 
process, while the other one satisfies neither of these properties. We 
prove that both are semimartingales and provide semimartingale de- 
compositions. 

Keywords: Levy motion, multistable process, semimartingale. 

*Universite de Rennes 2 - Haute Bretagne, Equipe de Statistique Irmar, UMR CNRS 6625, 
Place du Recteur Henri Le Moal, CS 24307, 35043 RENNES Cedex, France 
**Regularity Team, Inria, Ecole Centrale Paris - Grande Voie des Vignes, 92295 Chatenay- 
Malabry Cedex, France 

ronan.leguevel@univ-rennes2.fr, jacques.levy-vehel@inria.fr, limngliu.422@gmail.coni 

1 Background on multistable Levy motions 

The class of multistable processes was introduced in and further developed 
in [SI El [71 El Ho]. These processes extend the well-known stable processes 
(see, e.g. by letting the stability index a evolve in "time". Such models 
are useful in various applications where the data display jumps with varying 
intensity, such as financial records, EEG or natural terrains. Three paths 
have been explored so far to define multistable processes: the first one uses 
a field of stable processes X{t, a), and obtains a multistable process by con- 
sidering a "diagonal" Y{t) = X{t,a{t)) on this field. This is the approach 
of 12 E]. In j5], multistable processes are obtained from moving average 
processes. Finally, in [Hllin], multistable measures are used for this purpose. 

In this work, we shall be interested only in multistable Levy motions, 
which are the simplest examples of multistable processes. These multistable 



extensions of a— stable Levy motion were constructed in [H [6] using re- 
spectively Poisson and Fergusson-Klass-LePage representations, while similar 
processes were defined in ^ ilOj by their characteristic functions. In order to 
give precise definitions, let us set the following notations, which will be used 
throughout the paper: 

• a : M — > [c, d] C (0, 2) is a continuous function. 

• a : R — i- M is a function. 

• n is a Poisson process on [0, T] x M, T > 0, with mean measure the 
Lebesgue measure C 

• (ri)i>i is a sequence of arrival times of a Poisson process with unit 
arrival time. 

• {Vi)i>i is a sequence of i.i.d. random variables with uniform distribution 
on [0,T]. 

• (7i)j>i is a sequence of i.i.d. random variables with distribution P(7i = 
1) =>(7, = -1) = 1/2. 

The three sequences (rj)j>i, (Vi)i>i, and (7i)i>i are assumed to be indepen- 
dent. 

The following process, that we shall call field-based multistable Levy mo- 
tion, is considered in [4j: 

LMt) = a{t)C'J;l'^ 1m(X)Y<-V-W> (tG[0,T]), (1) 

{X,Y)Gn 

where y<-V"W> ;= sign(y)|r and 

Cu = (^j^ a;^"sin(x)dx^ . (2) 

Note that, when a{t) equals the constant a for all t, LpB is simply the 
Poisson representation of a— stable Levy motion. If one considers instead the 
Fergusson-Klass-LePage representation of Levy motion, then it is natural to 
define L fb SiS follows: 



+00 



LMt) = 5^7,rrV-Wij^^^,(y^) (t e [o,t]). (3) 



i=l 



This is the approach of [6], where it is proven that ([T]) and ([2D indeed 
define the same process, and that the joint characteristic function of Lps 
equals: 
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E I e^- ^ ^ ) = exp I -2 ^^^^ sin' \y2eMtj)^^^lioMi^) | dy dx 



(4) 

where m G N, (^i, ...,e^)e M'", (ti, . . . , t^) e M™. 

In [HI [To] another path is followed to define a multistable extension of Levy 
motion. Considering the characteristic function of a— stable Levy motion L: 

E(e'^^W) =exp(-t|e|"), 

one defines the process Ljj by its joint characteristic function as follows: 

ds I . (5) 



E I exp I z ^ e^Ljj{t^) j j = exp - y 



[o,t 



It is clear from this definition that Ljj has independent increments (this is 
why we call this version the independent increments multistable Levy motion). 
This is a strong difference with Lps, which is not even a Markov process. 
The construction based on multistable measures thus offers the advantage of 
retaining some properties of the classical Levy motion: obviously, increments 
stationarity cannot hold anymore, but we still deal with an additive process. 
This is not the case for the field-based Levy motion. However, a drawback 
of (|5]) is that, while LpB coincides, at each fixed time, with a Levy motion, 
this is not true of Ljj. 

The aim of this work is to elucidate the links between LpB and Ljj 
(section [2]), to prove that both processes are semimartingales and to give 
semimartingales decompositions of these processes (section [3]). 



2 Series representations of independent incre- 
ments multistable Levy Motion 

While three characterizations of LpB are known, namely the Poisson and 
Fergusson-Klass-LePage representations and its characteristic function, only 
the latter is available for L//. In this section, we provide the series represen- 
tation of independent increments multistable Levy motion. 

2.1 Poisson series representation 

Proposition 1. Ljj admits the following representation in law fort G [0,T].' 

Lnit)=ait) Cl{'ifh[oMX)Y<-'/'^^''>. (6) 

(x,y)Gn 
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Proof. To simplify, we take a{t) = 1 in this proof. Call L the process on the 
right hand side of (E]). Proposition 4.2 in entails that the series converges, 
and that its marginal characteristic function reads: 



E(e^^^W) = exp (-2 J J sin^ (^<{.t^l[o,](^)|l/r^/"(^)) dsdy^ 



(an integral sign without bounds means integration over the whole domain). 
We first prove that this quantity is equal to exp ^— \9\°'^'^'> ds^ . One com- 
putes: 

1 2sin^ (^<it^l[o,t](3)|y|-^/"(^)) dy = | (l - cos(<{J^)l[o,](.)br^/"(^))) dy 



4s) 



where we have used the change of variables z = 60^^^^'^^^'^ liQ^t]{s)y and 
the fact that 

f (1 - cos(2:)) 



\zH')+^ C, 



a{s) 



Thus, 

E(e^^^W) = exp(^-| 1 2 sin^ Q<(;(^)l[o,](.)y-^/"(^)) rf.rfy) 
= exp f- / |0|"(^)c/s 



The process defined by the series on the right hand side of (jnD clearly has 
independent increments. This fact and the above computation ensures that 
it coincides in law with Ljj. 

□ 

2.2 Fergusson-Klass-Le Page series representation 

To simplify the notation, we take T = 1 in this section. 
Proposition 2. Ljj admits the following representation in law: 

oo 

Ln{t) = ait) J2 C^«{^r^7.r-^/"(^'^l(v.<.). (7) 

i=l 
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Proof. We set again a{t) = 1. Let us first prove that the marginal character- 
istic function of the right hand side F of ([7]) is 

E(e'^^(*)) = exp . (8) 

Set 

where mx = {'y,U). Using the Marking Theorem (see [31 p. 55]), since 

P(7 = 1) = P(7 = -1) = 1/2, one computes: 



= exp J{1- e-'<(")"'^^-'^"*"'i("<*))(ixP(x, dm) 
= exp (^-i y /~ ((^ - e.ex-v.wi(„<,, j ^ _ ^^^^^ 



exp I — / / I 1 I dxdu 




'0 ^0 

fOO POD 




exp I - / / (1 - cos(eci(°|")l(„<t)X-i/"(")))dxdnj (9) 
= exp(^-y (10) 

where the passage from (jH]) to (fTOl) follows along the same lines as in the 
proof of Proposition [1] Now, 



E (e-^(^W-^W)) = exp(^y y (1 - e-^<"t'^^"''^'"'^(-"^'))rfxP(x, rfm) 

= exp (^-ly y~ ((l - e*^--^'°'"'i(=<"<*)) + (l - e-^-~^'"'"'i(-"<*))) da;rfn^ 



exp I — / / I 1 I dxdu 




'0 ^0 

f OO POO 



exp(^-y^ y^ (l-cos(^Ci(::|")l(,<„<,)X-i/"(")))rfxrfw 
expf- / |l(,,i](M)^r(")rfM 





ft 

exp ( - / 



This entails independence of the increments and ends the proof. □ 
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3 Multistable Levy motions are semimartingales 

3.1 Case of independent increments multistable Levy 
motion 

Proposition 3. Ln is a semimartingale. 

Proof. We use Theorem 4.14 in chapter II of [9]. Indeed, Ljj is an in- 
dependent increments process, and it suffices to verify that the function 
t E(e*^^'-*'') has finite variations on finite intervals for each 9. This is 
obvious in view of ([S]). □ 

It is easy to check that the system of generating triplets (in the notation 
of [T^) of Ljj is (0,z/, 0), with u{dx,dz) = \z\~°'^^^~^dzdx. If a ranges in 
[c,d] C (0,1), then Lu is a finite variation process, while if a ranges in 
[c,d] C (1,2), then Izlu^dXjdz) = {zl^^^^^dxdz < oo, which 

entails that Ljj is a martingale. The proof is a simple adaptation of the 
one for stable processes and is left to the reader. In general, the following 
decomposition holds: 

Proposition 4. 

Ln{t)=A{t)+M{t), 

where 

A{t) = a{t) C7i(xri[o,](X)r<-V"W> 
(x,y)en,|y|<i 

is a finite variation process and 

M{t) = a{t) <xr^l[o,](X)y<-V"W> 
(x,y)en,|y|>i 

is a martingale. 

Proof. That A has finite variation is a direct consequence of the fact that it 

is almost surely composed of a finite number of (jump) terms. 

To prove that M is a martingale, it is sufficient to show that it is in Li{Q) for 

all t. Note that the jumps of M are bounded by = sup^gjo t] ^i^) ^^Pbe[c d] ^1^^ • 
As a consequence, with obvious notations, 

pt poo pt pK pt pK 

/ / \z\vM{dx.,dz) < I I \z\i'{dx,dz) = 2 / / {z]'"'^^^ dx dz < oo. 
Jo Ji Jo Ji Jo Jl 

This entails that, for all t, M{t) has finite mean (see, e.g. [T?, Theorem 25.3]. 

□ 

Remarks 
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• A corresponding decomposition holds of course for the Fergusson-Klass- 
LePage representation, i.e.: 

Lii{t) = A'{t) + M'{t), 

where 

oo 

A{t)-a{t) C^fjj^-^ 7iii Hu,<t) 

j=i,ri<i 

has finite variation and 

oo 

j=i,ri>i 

is a martingale. 

• As is well-known, the decomposition above is not unique. Another 
decomposition of interest is the following: 

Ljj{t) = Mi{t) + A^{t) 

with 

oo 

i=l,a{U,)>\ 

is a martingale and 

oo 

^=1MU^)<\ 

is an adapted process with finite variations. 

3.2 Case of field-based multistable Levy motion 

LpB is not an independent increments process, and, contrarily to what its 
definition might suggest at first sight, it is not a pure jump process (see 
below for a more precise statement). Thus it is not immediately obvious 
that it is indeed a semimartingale. We shall use the characterization of 
semimartingales as "good integrators" (see, e.g. [21 [13]) to prove this fact. 
More precisely, fix t > and consider simple predictable processes ^ of the 
form 

n 

^(m) = ^o1{o}(m) + 

k=l 
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where < si < ti < S2 < t2 < ■■■ < Sn < tn = t, G J-'s^. and < 1 a.s. 
for all < k < n. The integral of ^ with respect to a process Y is: 

n 

fc=i 

It is well known that F is a semimartingale if and only if the family {/y(^), |^| < 
1,.^ is a simple predictable process} is bounded in probability. 

In fact, we shall establish the semimartingale property for a more general 
class of processes defined through Ferguson-Klass-LePage representations as 
follows: 

oo 

X{t) = a(t)C^g^*^ 5^7.r-^/"(*V(t, V^, (11) 

i=l 

where the function / is such that, for all t, |/(t, x)!"*^*-* dx < oo, and give 
conditions on the kernel / ensuring that X is a semimartingale. 

Theorem 5. Let X be defined by fllip . with a a function. Assume that X 
is a cddldg adapted process, and that there exists a positive constant such 
that for all {t,x) G M x [0,T], \f{t,x)\ < Loo- Assume in addition that, for 
all X E E, the function u t— ?■ f{u,x) has finite variation on finite intervals, 
with total variation Vf{.,x) verifying 

I \Vf{.,x)fdx<+oo (12) 
for some p G [d, 2). Then X is a semimartingale. 

Proof. We will show that itWt ^ semimartingale. Both processes t 

7ir^ and t h-> f{t,Vi) have finite variation, thus the same holds for 

OO 

t ^ ^ir~^^''^'^f{t,Vi). We shall prove that Y{t) = E 7*rr^^°^*V(^, ^i) is a 

i=2 

semimartingale. 



8 



One computes 

k=l 
n +00 



k=l i=2 

n +00 



fe=l i=2 
n +00 



fe=l i=2 



Let K > 0. We need to control P{\Iy{C)\> K). 

For the term A, we shall use the mean value theorem: there exists a sequence 
of random variables w'l G {sk,tk) such that 



(*.-.)^(iogr,)r-""«). 



One has 



P 



P 



n +00 



k=l 1=2 

(+00 n , I ( k\ 

1=2 k=l ^ \ I J 



< — E 
- KP 



t-00 n 



E E ( ^ ^^'> 



i=2 k=l 



(iogr,)r, 



As in the proof of Proposition 4.9 of |7], we use Theorem 2 of [Ij (note that 
random variables are independent of the {'Jj)j)- Since p<2, 
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p 



n +00 



fc=l i=2 



< —Ye 



n y I ( k\ 



+00 



< 



i=2 
+00 

i=2 



7i 



k=l 



> K 



PI 



.fc=l 



Kc^ 



p +00 



i=2 



where the infinite sum in the last fine above converges since p > d. Let us 
now consider the second term 5 of /y(^) : 



P 



n +00 



< — E 
- KP 



i=2 
+00 n 



> K 



i=2 k=l 



We use again Theorem 2 of [Ij : 

n +00 



< 



p 

2 
2 



k=l i=2 
+00 



> K 



i=2 

+00 

< —Ye 

KP ^ 

i=2 



< 



KP 
2 



+00 



i=2 

+00 



ifitk,Vi)-fisk,Vd) 

k=l 
n 

n 



-L ~r J- 7 



i=2 



+00 



KP 
2 



i=2 



-1/c ^ p-l/d 



+00 

i=2 



1/d 
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where, in the last hne above, the first expectation is finite by assumption 
(fT^ and the infinite sum converges converges because p > d. 

We have thus shown that P{\Iy{0\ > K) < ^ for some constant C, as 
required. □ 

Corollary 6. Assume a is a function. Then the field-based multistable 
Levy motion is a semimartingale. 

The result above does not give a semimartingale decomposition of Lps- 
The following theorem does so. In addition, it elucidates the links between 
the two multistable Levy motions: 

Theorem 7. Assume a is a function. Almost surely, for all t, 

LpB{t) = A{t) + Ln{t), (13) 
where A is the finite variation process defined by: 



A{t)= / ^7. ^ '\ Hs)lio,smds. (14) 

"^0 i=l 

From a heuristic point of view, this result states that, at all jumps points, 
both Levy multistable motions vary by the same amount. "In-between" 
jumps, however, Ljj does not change (it is a pure jump process), while LpB 
moves in a continuous fashion (and thus is not a pure jump process). 

Proof. For the sake of notational simplicity, let: 

9^{t) = c]!;^'T-''-^' = h{t)T-''-'\ 
dm = m = mr;'^"'' + ^«rr^/"(*) inr., 

e,{t) = 6'(t)r^/"W + MM^)ri/"W(lnz). 

a'^{t) 

We wish to prove that (again, we set a{t) = 1): 

+ 00 „4 +00 +0O 

i=l 1=1 i=l 

Let 

N 



Dn{s) = E7i(ii(s)l[o,s[(^i), 



i=l 
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and 

+00 

i=l 

Finally, write: 

N 

En(s) = J^7iei(s)l[o,.[(V^i), 
1=1 

and 

+00 

-^(*) = X^7iei(s)l[o,5[(l^i). 

i=l 

The function s 1— )■ 7i(ii(s)l[o,s[(Vi) is continuous and difFerentiable on (V^, 1), 
and vanishes on [0, l^]. 

Choose AeR such that Vt e (0, 1), \b'{t)\ < A , \^^^^^\ < A and < 
A. 

Case d < 1: 

Fix 1 > > and io e N such that, for i > io, ^^JJ^ < r^. Then, for all 
s e [0, 1] and p > po, 

g 

i=p 

< 



< 
< 



The uniform Cauchy criterion entails that Er^{s) converges uniformly on [0, 1] 
to E{s). 
Case d > 1: 
Set 

23+l_i 
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and b{j) = j^^^^^^0^V2^. Consider 

m - 



sup ms)\<j^^^±il±^V2. 



SG[0,1] 



sup \Ei{s)\<b{j) 

5G[0,1] 



Fix d > d > 1 and io e N such that, for i > io, ^^^l^P/'"' < Let 



5 > 2 + ^ ~ 2- Define the step function aj{s) — '^fiil*)- 
For all s e [0,1], 



k 



\a 



< 



fe=0 

A 



and thus 



Set 



sup \aj{s)-a{s)\ < j-^. 
se[o,i] i^^ J 



i=23 



Then, 



23+1-1 



2J+l_l 



i=2J 



The finite increments theorem entails that there exists a^{s) € [aj{s),a{s) 
(or in [a{s),aj{s)]) such that 



2-'+i-l 



In 2 



i=23 



a{s)a'{s 



(a[s) — aj(s 



'i'+^-l /, •\2 



=21 



13 



and thus 



A 

< 



2 



< 



Choose Jo e N such that, for all j > Jo, 

sup \Ei{s)-E^s)\ < ^. 

Then, for j > Jq, 

p(e(j)) = P( sup \El{s)\>h{j) 
^ ' \se[o,i] 

< p(M + sup \Ei^{s)\>h{3) 

= P( sup \Ei^{s)\>^-^^^^ 
\se[o,i] 

Denote the set of possible values of OLj. Then Cardi^A^^ = [2^^]. For each 
random drawing of the {Vi)i, E^.{s) is composed of a sum of n terms of the 

form where ckq £ . There are 2^ possible values for n, and [2^''] 
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possible values for Oq. One thus has the following estimates: 

b{j) 



f(e{j)) < P sup \Ei {s)\> 



2^ ,\r^ ,-i/«n, b{j) 



< P ^i,,...,i,M2^,2^^^-ii{Vh <Vi,< ... < Vijn { sup sup I 5^7/.^r'^"°l > ^} 

k 

< E p ( {^/, < < ... < } n { sup sS'p I ^ 7/./r'^"° I > ^} 



/l,...,i«G[2J,2J+i-ll 



< Yl ^{vi,<vi,< ... < vi^j p sup sup > 

ii,...,i2,e[2^2^+i-l] \ao(^A^ k-1 .^^ 

/l,...,i2jG[2J,2J+i-lI ^ \ i=l 



where Stute's lemma |14) was used in the end. As a consequence, 

P ^liminf = 1. 

Let now p G N, and set jp = [j^]- Define 

2JP+1-1 

K^'is)= E 7.e.(s)l[o,.[(V^.), 

i=m 

and 

As above, 

sg[0,l] 



sup - a(s)| < 



15 



and 



^2 2^^+'-l 



\E'^'P('i) — E'^'P('i)\ < \ 



< 



< 



Fix po such that for all p > po, 



mel2Jp,2*+i-l] se[o,i] 

and consider 



sup sup \E^'P{s) - E^j^{s)\ < 



Eip) = <! sup sup \E:^'P{s)\<b{jp) 

me|2Jp,2Jp+i-ll se[0,l] 



For p > Po, 

KJp) 



P 



(e{p)) < p sup sup |^:::;^(.)|> _ 

^ ^ \m6[2^f,2^P+i-l]se[0,l] ^ 



< 

mel2-'P,2JP+i-l] 



me|2^P,2Jp+i-l] 

As a consequence, 

P ^liminf = 1. 

Reasoning along the same lines and setting 

m 

K'Ks) = E 7ie.(«)l[o,.[(^0, 

i=2^9 
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and 



Eiq)^{ sup sup |£;r(«)|<6(i,)^, 

mel2^9 ,2^9+1-1] s6[0,l] 



one gets 



P l^imM E{q)j = 1. 
Finally, let p, g e N, with q > p: for all s e [0, 1], 

lie-i{s)lio,si{Vi) + ^ ^ 7iei(s)l[o,s[(^i) + 5Z 7iei(s)l[o,5[(^0 

< b{j,)+ b{j)+b{j,) 

+ 00 

j=jp 

Again, the uniform Cauchy criterion entails that En{s) converges uniformly 
on [0, 1] to E{s). 

Convergence of the difference term: 

Let q>p and denote Rp,q{s) = J2 li{di{s) - ei(s))l[o,s[(V'i)- 

i=p 

Fix po such that, for all p > po, 



One computes: 



Rp,q{s) < -ei(s)|li. 



< |i-v«(*) - r-^/"(*)|i,<r,<^ + inz - rr^/"W inr,|i,<r. 



<2 



9 



^ zV<=«' Inz S 



< 



(l + lni) 



— - 1 
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The series ^ '-^^j,"^^ | ^4 — 1 1 converges almost surely. The uniform Cauchy 

i 

criterion thus applies to the effect that Dn{s) — En{s) converges uniformly. 
As a consequence, Dn{s) converges uniformly to D{s). 
Thus, almost surely. 



Besides 



pt pt pt +°° 

lim / DN{s)ds = / D{s)ds = / 7i5f-(s)l[o,s[(V"i)ds 
^+~Jo Jo Jo ~[ 



DN{s)ds = g'i{s)l[o,s[{Vi)ds 
i=i Jo 



N 



= Y.^^{g^{t)-g^ml[0,tm 
1=1 

N N 

= J2^,g,{t)^odV^ ~J2^igi{Vi}lio,t]{V^ 

i=l i=l 

which entails finally that 

pt +00 +00 

lim / Div(s)cis = V7i^i(t)l[o,t](l^i) - y'7i^i(V^»)l[o,t](^i). 

1 = 1 1=1 



This is (USD. 

That A has finite variation follows from the fact that it is an absolutely 
continuous function. □ 
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